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Charge carrier mobility in disordered organic materials is being actively studied, motivated by 
several applications such as organic light emitting diodes and organic field-effect transistors. It is 
known that the mobility in disordered organic materials depends on the chemical potential which in 
turn depends on the carrier concentration. However, the functional dependence of chemical potential 
on the carrier concentration is not known. In this study, we focus on the chemical potential in 
organic materials with Gaussian disorder. We identify three cases of non-degenerate, degenerate 
and saturated regimes. In each regime we calculate analytically the chemical potential as a function 
of the carrier concentration and the energetic disorder from the first principles. 
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Charge transport in disordered organic materials is ex- 
ploited in a wide range of devices, including organic light 
emitting diodes (OLEDs)ji. organic field-effect transistors 
(OFETs)^ photoreceptors^ and photovoltaic cells^ In 
disordered organic materials, the carrier mobility is due 
to thermally assisted tunneling "hopping" between lo- 
calized molecular states££ It is know that the carrier 
mobility depends on the temperature, energetic disorder, 
and carrier concentration.^ For efficient device model- 
ing it is useful to have a compact analytical expression 
for mobility. An attempt was made by Coehoorn et. alS- 
where they provided an analytical expression for mobil- 
ity in organic materials with Gaussian energetic disorder. 
However, the authors noted that the expression for mobil- 
ity is not suitable for practical numerical device modeling 
because no analytical expression for the chemical poten- 
tial as a function of carrier concentration and energetic 
disorder is available. The main purpose of this paper is 
to show that an analytical expression for the chemical 
potential in organic materials with Gaussian energetic 
disorder can be obtained from the first principles, with 
no free parameters. We show that the derived expression 
for the chemical potential is fairly accurate and the error 
involved is well below the thermal energy. 

Within the Gaussian disorder model, it is assumed that 
the density of states (DOS) is given by: 

where a is the standard deviation of the DOS and is a 
measure of the energetic disorder while iV t is the total 
number of hopping sites per unit volume. For a given 
charge carrier concentration, p, the chemical potential, 
/i, is related to the density in the following way: 
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where fee is the Boltzmann's constant and T is the tem- 
perature. 

In general, for given p and c, Eq. ((2]) is solved itera- 
tively to obtain the chemical potential. However, in the 
limit of vanishing carrier concentration, an analytical ex- 
pression for the chemical potential can be easily obtained 
because in the limit of small p carriers can be considered 
essentially independent of each other. In this case Eq. ([2]) 
can be solved for /j, by replacing the Fermi-Dirac (FD) 
statistics with the Boltzmann statistics^ In this limit, 
referred to as the Boltzmann approximation (BA), the 
density of occupied states (the product of DOS and the 
FD distribution function, DOOS) is to a very good ap- 
proximation given by a Gaussian, centered at the energy 
value E = —tj 2 /(k-B,T) .— The corresponding expression 
for /i is given by: 

/J = -777-^ + k B T In p. (3) 
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The BA limit is applicable to OLEDs, where under typi- 
cal operating conditions, the concentration is 10~ 4 -1CP 5 
carriers per hopping site<£ On the other hand in OFETs^ 
application of high gate voltage can lead to a concen- 
tration of 0.01-0.1 carriers per hopping site. At these 
carrier concentrations, interaction between carriers be- 
comes significant implying that Pauli's exclusion princi- 
ple must be taken into account. In this high concentra- 
tion regime, the Boltzmann approximation is no longer 
valid and Fermi-Dirac distribution must be used. We 
show that even in this regime, an analytical solution for 
\x can be obtained in a simple and intuitive manner. In 
the following text we identify three regimes referred to 
as the non-degenerate, degenerate, and saturated regime. 
In each regime we calculate analytically the chemical po- 
tential as function of the carrier concentration. 

Our objective is to solve Eq. ([2]) for the chemical po- 
tential, /i, for given values of carrier concentration, p, en- 
ergetic disorder er, and temperature T. We evaluate the 
integral in Eq. @ using a saddle point approximation. 
We assume that the saddle point, E*, of the integrand 
is known. Assuming this one gets two coupled algebraic 
equations for the chemical potential, [i, and the saddle 
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point of the integrand instead of one integral equation for 
the chemical potential^ The first algebraic equation is 
obtained by imposing the condition that E„ is the saddle 
point of the integrand in Eq. © . On doing so we get the 
following equation for the chemical potential 



fx = E* + fc B Tln ( — 



E*knT 



- 1 



(4) 



Saddle point approximation implies that the integrand 
(cquivalently the DOOS) is a Gaussian centered at E % 
with a standard deviation ov Namely, Eq. ^ is ap- 
proximated by 
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where the expression for the DOOS is 

N t e~^ (E*k B T + <j 2 
h(E) = )= '- exp 



2na 3 



{E-E* 



2a 



and 



fc R T(T 4 



a 2 (k B T- E*) - E 2 k B T' 



(6) 



(7) 



On evaluating Eq. (0 we obtain the second algebraic 
equation: 
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The two coupled equations Eq. (Q| and Eq. (5) can be 
solved to obtain fi and E*. We now solve Eq. Q and 
Eq. ([8]) self-consistently in three limiting cases: 

a. Non-degenerate regime: This regime corresponds 
to the case of vanishing carrier concentration. Eq. ([8]) 
implies that in this regime the saddle point occurs at 
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The dependence of E* on the carrier concentration p is 
obtained as 
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Self-consistency in this regime requires that carrier con- 
centration is low enough such that the condition ([9]) 
holds. This implies that in this regime p <C Pi, where 
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The chemical potential in this regime is well known and 
is given by 
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Figure 1: Chemical potential fi vs. the normalized carrier 
concentration p for different values of a/ksT (numbers on 
the plot). Full lines and dashed lines correspond to the non- 
degenerate, Eq. (|12|l . and degenerate, Eq. (15) , regime, re- 
spectively. Symbols correspond to the numerically calculated 
chemical potential. 



The non-degenerate regime corresponds to the BA limit. 
In this regime an increase in carrier concentration leads 
to an increase of the maximum of DOOS with negligible 
shift in E*& 

b. Degenerate regime: For p 3> p\ the saddle point, 
E*, shifts significantly from —a 2 /k B T and the BA limit is 
no longer valid. In this regime the carrier concentration 
is such that saddle point satisfies —<r 2 /k B T <C £7* <C 




Figure 2: Comparison of analytical result, Eq. (|19[) (lines), 
and the numerical calculation (squares) . Full lines correspond 
to the degenerate regime. The horizontal dotted line corre- 
sponds to the limit of the non-degenerate regime. 
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This analytical expression for the chemical potential in 
the degenerate regime is the main result of this paper. 
Since our focus is on the systems with a > k B T, the 
second crossover value, P2 is of the order of one and de- 
generate regime is valid until almost the full saturation 
of all the available states. For the sake of completeness 
we discuss below this saturated regime. 

c. Saturated regime: In this regime the carrier con- 
centration, p, approaches the limiting value of pi ~ 1 
where almost all the states are occupied and the DOOS 
is the same as the density of states. In this limit, the sad- 
dle point satisfies > E* ^> — k B T, so Eq. (|SJ) implies 



E* = -o^/21n- 
P 
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Self-consistency in this regime requires p > p2- The 
chemical potential in this regime is given by 



(jlll) . with increasing a, the crossover from the non- 
degenerate to the degenerate regime occurs at lower con- 
centrations. The crossover from the degenerate to the 
saturated regime is not shown in the figure. The satu- 
rated regime, irrelevant for organic devices, is a narrow 
regime with divergence of chemical potential at p = 1. 
It is clear from the figure that the analytically obtained 
chemical potential and the crossover density, Eq. (jlll) . 
are in good agreement with the numerical calculations. 

In addition to the mobility of carriers, to model charge 
transport one needs to calculate the diffusion coefficient. 
Assuming equilibrium conditions, the diffusion coefficient 
is determined from the the generalized Einstein relation^ 
which is given as 
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This regime is irrelevant for organic devices because even 
in OFETs, where carrier concentration can be tuned 
by applying gate voltage, the concentration is typically 
much smaller than l£ 

In Fig. 1, we plot the chemical potential as obtained 
above for different values of a. As predicted by Eq. 



where 7 is ratio of diffusion coefficient to the mobility 
of carriers and q is the elementary charge. Roichamn 
and Tessler— showed that in a Gaussian DOS, 7 de- 
viates significantly from the low density limit value of 
k B T/q. They obtained 7 as an implicit function of the 
chemical potential, /i. Our expression for the chemical 
potential allows us to obtain an analytical expression for 
the generalized Einstein relation in each of the regimes 
as a function of p, a and T: 
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In Fig. 2 we compare our analytical expression for 7, Eq. 
(fT9|) . with the numerical results (notice that the inverse 
of 7 is shown). As expected, a fairly good agreement is 
obtained far from the crossover concentration, p\. 

The results above indicate that the analytical expres- 
sions for the chemical potential are fairly accurate. Nev- 
ertheless, it is useful to estimate analytically the error 
in the derived formulas. The main source of the error is 
the saddle point approximation, Eq. (|6|), in evaluating 
the integral in Eq. @. Under the saddle point approxi- 
mation, the integrand (or equivalently the DOOS) of the 
integral in Eq. @ is a Gaussian centered at E„. The in- 
tegrand is well approximated by a Gaussian only in the 
neighborhood of E* (within few standard deviations, cr*). 
Away from E tl the deviation of the integrand from the 



Gaussian gives rise to the error in evaluation of Eq. © . 
The relative error in saddle point approximation of the 
integral can be written as 

00 
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For energies well below the saddle point, E -^i E* — 
3c*, the FD distribution can be approximated as unity 
whereas well above the saddle point, E ^> E* + 3er*, 
it can be approximated by the Boltzmann distribution. 
Therefore the error can be estimated as 
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The integrals appearing in the above expression can be 
simplified using error function. For given values of car- 
rier concentration, p, and energetic disorder, a, the other 
parameters E t , er*, and /i can be calculated as shown 
above. In Fig. 3 we show that for relevant range of a 
and p values, the value of 5 does not exceed 20%. This 
implies that the error in the chemical potential is well 
below than the thermal energy, k^T. 

Although we derived the expressions for the chemical 
potential only in the respective regimes, the algebraic 
equations Eq. and Eq. © can be easily solved nu- 
merically to any desired degree of accuracy for any given 
values of p and a. 

To conclude, in this paper we derive from the first prin- 



ciples an analytical expression for the chemical potential 
in organic materials with Gaussian disorder. In the most 
relevant, degenerate regime our result, Eq. (|T5|) . is new 
and, as we demonstrate, is fairly accurate. We show that 
over the relevant range of carrier concentration and en- 
ergetic disorder, the error in calculation is well below the 
thermal energy. In addition, the existing iterative numer- 
ical techniques to calculate chemical potential can use our 
derived expressions as fairly accurate starting point. This 
can lead to much more efficient algorithms to calculate 
chemical potential and therefore the mobility in organic 
materials with Gaussian disorder. 

We thank M. Depken and P. A. Bobbert for helpful 
discussions. 
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